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ABSTRACT 

We investigate the effects of radiation pressure from stars on the survival of 
the star-forming giant clumps in high-redshift massive disc galaxies, during 
the most active phase of galaxy formation. The clumps, typically of mass ^ 
10^ — lO^Af© and radius ^ 0.5 — 1 kpc, are formed in the turbulent gas-rich discs 
by violent gravitational instability and then migrate into a central bulge in 
^ 10 dynamical times. We show that the survival or disruption of these clumps 
under the influence of stellar feedback depends critically on the rate at which 
they form stars. If they convert a few percent of their gas mass to stars per 
free-fall time, as observed for all local star-forming systems and implied by the 
Kennicutt-Schmidt law, they cannot be disrupted. Only if clumps convert most 
of their mass to stars in a few free-fall times can feedback produce significant 
gas expulsion. We consider whether such rapid star formation is likely in high- 
rcdshift giant clumps. 

Key words: galaxies: formation — galaxies: ISM — galaxies: star clusters: 
general — galaxies: star formation — ISM: clouds — stars: formation 



1 INTRODUCTION 

A significant fraction of the massive galaxies, ~ 
10^^ Mq in baryons, during the period from z = 
1.5 — 3 when star formatio n is at its peak and most 
stellar mass i s asse mbled (|Hopkins fc BeacomI l2006l : 
iMagnelh etlj] |2009| ). form stars at high star forma- 
tion rates (SFR) of ~ lOOM0yr"\ Many of these are 
turbulent, gas-rich, extended rotating discs in which 
much of the star formation takes place in a few gian t 
dumps JCowie et al.l [l995l: Ivan den Bergh et all 1 19961 : 



lElmegreen et al.ll2004al . i2005l . l2007l ). Based on this mor- 
phology, they were first termed "chain" or "clump- 
cluster" galaxies. In a typical galaxy of this type, 
10 — 40% of the UV rest-frame light is emitted from a 
few clumps of characteris tic size ~ 1 kpc that form stars 
at tens of Mq yr~^ each (lElmegreen et a"l]|2004bl . l2005l : 



clumps, and their survival subject to stellar feedback, 
are central to our understanding of galaxy formation. 

Kinematically, most of the galaxies that host the 
clumps are thick rotating discs, with high velocity dis- 
persions of cr = 20 — 80kms^^ (one dimensional), 
compared to cr ~ lOkms^^ is present-day discs; they 
have rotation to dispersion ratios of V/a ~ 1 — 7 
(|Cresci et al.ll2009l ). Estimates of the total gas fraction 
in star-forming galaxies, based on CO measurements, 
range from 0.2 to . 8, with an average of ~ 0.4 — 0.6 
(|Tacconi et al.ll2008l : lOaddi et all [20081 . 120091 ). system- 
atically higher than the typical gas fraction in today's 
discs. These properties are generally incompatible with 



iForster Schreiber et al.ll2006l : lGenzel et al.boO^ '). These 
clumps are much more massive than the star-forming 
complexes in local galaxies. The star formation in these 
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Instead, gia nt clumps f o rm thr ough a scenario, 
summarized by iDekel et al.l (|2009bl ). in which mas- 
sive galaxies at z ~ 2 — 3 are fed by a few nar- 
row and partly clumpy streams of cold gas (~ 10*K) 
that flow along the dark matter filaments of the cos- 
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mic web (e.g. iHahn at all I2OO1I ) and penetrate deep 
into the centres of the massive dark matter haloes of 



~ lO^^Mp) jBirnboim & Dekel 20031: 


Keres et all 20051: 


Dekel & Birnboimll2006l: Ocvirk et al. 


20081: Dekel et al.l 


2009al). Indeed, the existence of these streams is an in- 



evitable prediction of the standard ACDM cosmology. 
They may pr oduce the structures we observe as Lyman- 
alpha blobs (iFurlanetto et al.ll2005l : iGoerdt et allbOOgl : 
iDiikstra fc Loebl |2009| ~ The angular momentum they 



carry leads the accreted material to form a disc of ra- 
dius Rd ~ lOkpc. The continuous intense input of gas 
at the level of ~ 100 Mq yr~^ maintains a high gas sur- 
face density E, which drives a violent gravitational in- 
stability with a Toomre Q parameter below unity, Q ~ 
aQ/{nG'S) < 1, where a is the one-dimensional velocity 
dispersion and fl i s the angular v elocity associated with 
the potential weU (|Toomrelll963 ). The disc instability is 
self-regulated at Q ^ 1 by the gravitational interactions 
in the perturbed disc, which keep the disc thick and 
with a high velocity dispersion. The disc forms strong, 
transient spiral features that fragment to produce 5-10 
bound clumps, that together comprise ~ 20% of the 
disc mass. The largest clumps have characteristic radii 
R ~ 7GE/f2^ ~ 1 kpc, and characteristic masses of a 
few percent of the disc mass, M ~ 10^ Mq . A spectrum 
of smaller clumps with somewhat lower masses forms as 
well. The clumps' large masses cause them to migrate to 
the centre of the disc on a short time scale of ~ 10 disc 
crossing times, where they rnerge in to a central bulge 
(lNoguchi|[l999l: llmmeh et al.ll2004allbl: iBournaud et al.l 



l2007l : ICeverino et al.ll2009l : lAgertz et al.ll2009t ). This vi- 
olent instability phase can last for more than a gigayear, 
during which the mass flow from the disc to the bulge is 
replenished by fresh accretion, keeping the mass within 
the disc radius divided quite evenly between disc, bulge 
and dark matter components. 

While this scenario of clump formation and migra- 
tion to build up bulges is appealing, it relies on the 
ability of clumps to survive for ~ 10 disc dynamical 
times, i.e. a few hundred million years, while the gas 
i n them turns into stars on a comparable timescale 
(iGeiizel et al.ll2008l : iDekel et al.]|2009bl : ICeverino et~aLl 
I2OO9I ). However, at the SFR of a few tens of solar 
masses per year in the clumps, it is possible that they 
might be disrupted by stellar feedback o n considerably 
shorter timescales. iMurrav et al] (|2009l ') argue for ex- 
actly this scenario. In their models, clumps disrupt af- 
ter ~ 1 dynamical time, during which they turn only 
~ 30% of their mass into stars. In this picture clumps 
would not survive long enough to migrate, and bulges 
would instead need to be built up by mergers. While 
this scenario seems difficult to reconcile with the es- 
timated ages of a few hundred Myr for the oldest stel- 
lar populations in some clump s (|Elmegreen et al ]|2009al : 
iForster Schreiber et al.l 120091 ), there is sufffcient uncer- 
tainty in both the observational estimates of clump ages 
and the theoretical modeling of clump evolution and dis- 
ruption to merit a re-investigation of the problem, which 
we provide in this paper. 



The outline of the paper is as follows. In ^ we 
derive the expected gas ejection fraction as a function 
of SFR efficiency. In ij3] we address the observational 
estimates of the SFR efficiency. In |4]we discuss some of 
the issues raised by our results and compare to previous 
work. We summarize our conclusions in 35] 



2 RADIATIVE FEEDBACK AND CLUMP 
SURVIVAL 

Consider a uniform-density giant gas clump of mass M, 
radius R, and surface density E = M/{-kR^). It forms 
stars at a rate Mt , and the stars formed within it have a 
combined luminosity L. (We defer discussing the effect 
of sub-clumping within giant clumps to § 13.31 since it 
does not change the qualitative result.) Characteristic 
numbers to keep in mind for the largest, best-observed 
clumps, found in galaxies with baryonic masses ~ 10^^ 
Mq, are M ~ 10^ Mq, R~1 kpc, and E ~ 0.1 g cm"\ 
Clumps in the ~ 10^" Mq galaxies, which are more com- 
mon, have masses ~ 10* Mq and sizes that are at or be- 
low the resolution limit of present observations. We will 
assume that they have surface densities comparable to 
those their larger cousins; they cannot be much smaller, 
since the mean column densities of the galactic disks as 
a whole is E ~ 0.05 g cm~^. We wish to evaluate the 
fraction e of clump mass that will be ejected by stellar 
feedback and the fraction £ = 1 — e that is transformed 
into stars, because this is the critical parameter that 
determines whether the clump will form a bound stel- 
lar syst em. Both N -body simulat ions and analytic mod- 
els (e.g.lHillslll980l:lKro upa 2001: iKroupa fc Boilvl[20o3 : 
iBaumgardt fc Kroupal 120071 ) indicate that if £ > 0.5, 
then most of the stellar mass will remain bound, while 
ii £ < 0.3 then no bound stellar system will be left. 
Small portions of the clump where £ was locally higher 
may form bound clusters, but these will be orders of 
magnitude smaller than the initial clump. 

Several common feedback mechanisms are not im- 
portant for giant clumps in the relevant mass range. 
First, supernova feedback is unlikely to be effective 
in ejecting mass from these gia nt clumps, due to 
cooling and leakage of ho t gas dDekel et al. | l2009bl : 
iKrumholz fc Matzneij l2009l : iMurrav et al.1 12OO9I ). Sec- 
ond, the pressure of warm (~ 10* K) ionized gas IS m- 
effective because the escape velocity from the clump is 
larger than the gas sound speed of ~ 10 km s~^. Third, 
protostellar outflows are unable to eject ma ss because 
they do not provide enough momentum (|Fall et al.l 
I2OIOI ). Instead, the dominant feedback mechanism is 
likely to be radiation pressure from newly-formed stars, 
which creates a radiation-dominated H 11 region. The 
expansion of such a region f ollows a similari t y solu tion 
(|Krumholz fc Matzneijl2009l '). and lFaU et al] ((2OI0I ') use 
this solution to show that all the remaining gas will be 
ejected once the fraction of gas mass transformed into 
stars reaches a value 
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where 



Ecrit — 



577/t.ap (L/M,) 
nacritGc 



(1) 



(2) 



77 is a constant of order unity that depends on the den- 
sity distribution within the clump, (L/M,) is the light to 
mass ratio of the stellar population, /trap represents the 
factor by which the radiation force is enhanced by trap- 
ping of re-radiated infrared light within the expanding 
shell, and Qcrit is a parameter of order unity that de- 
scribes the critical ve l ocity required to eject mass from 
the cloud. iFall et all (|2010l ) assume fiducial values of 
ri = 2/3 and /trap = acrit ~ 2 for the local star- forming 
regions, and we adopt the same values here since there 
is no obvious reason for them to be systematically dif- 
ferent in the case of giant clumps at redshift 2. We refer 
to £ as the final star fraction^ 

To determine the light to mass ratio, we must 
deal with the complication that the characteristic cross- 
ing time of a giant clump is rather long, ta = 
R/^GM/R = 1^M~^''^R\'^ Myr, where Mg = 
M/[l{f Mq) and Rx = R/(l kpc). Depending on the 
exact values of M and R, this can be either greater 
than or less than the main-sequence lifetime of mas- 
sive stars. Thus we can neither assume a single burst, 
so that all stars are coeval, nor continuous star forma- 
tion, so that the population is in equilibrium between 
new stars forming and old ones evolving off the main 
sequence. However, we can treat these scenarios as two 
limitin g cases, which must brac ket any real stellar popu- 
lation. iKrumholz fc TanI (|2007l ) point out that in stellar 



populations younger than 3 Myr, where no stars have 
left the main sequence yet, the light-to-mass ratio has 
a consta nt value, ^, for which we adopt ^ « 2200 erg 
s~^ g~^ (|Fall et al.ll2010l ). Once a stellar population is 
old enough to reach statistical equilibrium between star 
formation and star death, it instead has a nearly con- 
stant luminosity-to-star-formation ra te ratio, which 
we ta ke to be $ « 6.1 x 10^^ erg g~^ jKrumholz fc TanI 
120071 '). For any realistic stellar population whose light is 
dominated by young, massive stars as opposed to old 
ones, the luminosity is roughly equal to the smaller of 
these two limits, and for simplicity we simply take the 
light-to-mass ratio to be 



L 



A/, 



(3) 



We refer to the first case as the "young stars" limit and 
the second as the "old stars" limit, since they repre- 
sent the opposite extremes of stellar populations that 
have undergone no evolution and populations that are 



^ Note that this is sometimes referred to as star formation 
efficiency as well, but we avoid the term efficiency because it 
does not have a standard meaning, and different authors use 
it for different concepts. 



old enough to have reached equilibrium between star 
formation and stellar evolution. 

Substituting these two light-to-mass estimates into 
equation ([1} gives 



£ — max 



1 + 



-'crit,y 



1 



dcp 



where 



— ^ ^ 1.2 g cm 



, (4) 



(5) 



is the critical density evaluated in the young-stars limit 
and tdop = M/Mt is the depletion time, i.e. the time 
that would be required to convert all of the gas into 
stars. The numerical eval uation in equatio n (O is for 
the fiducial parameters of iFall et al.l l|2010t ). Note that 
this expression is a maximum rather than a minimum 
because £ is a decreasing function o f {L/M ,). 

Following iKrumholz fc McKed (120051 ). we define 
the dimensionless star-formation rate efficiency as the 
ratio between the free-fall time and the depletion tim^E 
namely 



(6) 



M/ta ■ 

IKrumholz fc TanI (|2007h show that eg ~ 0.01 across 
a very broad range of densities, size scales, and envi- 
ro nments. He re we define the free-fall time as ts = 
•y/37r/(32Gp), where p is the gas density. We discuss 
in § |3] below whether this value of es applies in high- 
z giant clumps. If we adopt it for now and make this 
substitution in equation (j4|, we find that 



£ = 



l-f 



1 + 



^crit.y 

E 
12 

e77 



-'crit.y 



,1 



(ttEM)!/-!* 

0.086 

(E_iA/9)i/4'ff'- 



(7) 



where E_i = E/(0.1 g cm ^), es^. 



Efi/lOO, and the 



numerical evalua tion uses the fiducial parameters from 
iFall et all (|2010h . We use Equation (jH]) to plot f as a 
function of es in Figure [1] 

What does this result imply for the survival of high- 
z giant clumps? We note that the second term in brack- 
ets in equation (|8]), corresponding to the case of a stel- 
lar population older than ~ 3 Myr, is generally the one 
that applies for giant clumps. This reflects the fact that 
the crossing time for our flducial values of the parame- 
ters is 15 Myr, which is signiflcantly larger than 3 Myr. 
The young stellar population limit applies only if the 



Note that the rate efficiency eg that we have defined here 
is distinct both from the star formation rate, which has units 
of A^Q yr~^ and is not normalized by the ratio of mass over 
free-fall time, and the final stellar mass fraction S, which is 
dimensionless but does not carry any information about the 
star formation rate. 
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< Old stars 1 I Young stors 




Figure 1. Final stellar fraction £^ as a function of star for- 
mation rate efficiency eg, computed from Equation |(7]|. We 
show results for giant clump surface densities S_i = 0.2, 
1.0, and 5.0 (green, red, and blue lines, as indicated), and for 
giant clump masses Mg = 0.03, 0.3, and 3.0 (dotted, solid, 
and dashed lines, as indicated). The arrows schematically in- 
dicate the region around efj = 0.1 where we change from the 
old stars limit (low eg) to the young stars limit (high eg). The 
gray region from £ = 0.3 — 0.5 indicates the rough boundary 
between stellar fractions £ < 0.3, for which no bound stellar 
system will remain, and £ > 0.5, for which a majority of the 
stellar mass will remain bound. 

star formation rate efficiency is much higher than is ob- 
served in any star-forming systems anywhere in the lo- 
cal universe, £fi,_2 > 10, or if the star-forming systems 
are significantly less massive or much more dense. This 
makes high-z giant clumps very different from Galac- 
tic star clusters or even super-star clusters, such as 
those found in local starburst galaxies, e.g. the An- 
tennae or M82. These have lower masses and higher 
surface densities, with cro ssing times ~ 0.1 Myr (e.g. 
iMcCradv fc GrahamllioOTi ). placing them firmly in the 
young-star limit. Indeed, the first term in the brackets 
is identical to that derived bv lFall et all l|20ld 'l for local 
star-forming clumps|f| 

We conclude that unless eft,_2 S> 1 for giant 
clumps, as opposed to ~ 1 for the local star- forming sys- 
tems, star-forming clumps with masses ~ 10* - 10^ Mq 
and surface densities ~ 0.1 g cm~^ cannot be disrupted 
by radiation pressure — so they end up converting most 
of their mass to stars. Thus, the expulsion fraction rel- 
evant for high-z giant clumps, as derived in the old-star 
limit, is 

e = 1 - f = 0.086EZY'^M^'^^'^es,-2 ■ (9) 
Since £ > 0.5, we expect the resulting stellar sys- 

^ To get numerical agreement in the coefficient of S_i, we 
must adjust fc by a factor of VO.4 to account for measuring 
the escape velocity at the surface instead of at the half-mass 
radius. 



tems to remain gravitationally bound. Clumps with 
E_i = 1 have £ < 0.5 and suffer significant disrup- 
tion only if their masses are below ~ 10® Mq , and they 
reach £ < 0.3 and undergo complete disruption only at 
masses < 2 x 10^ Mq. Thus the observed ~ 10*^ - 10^ 
Mq giant clumps in high— z galaxies should survive dis- 
ruption unless eff,_2 3> 1. Even if we are maximally 
conservative and assume that the smaller clumps have 
E_i = 0.5, i.e. that their surface densities do not exceed 
the mean surface densities of their host galaxies, our es- 
timated maximum mass for disruption only increases by 
a factor of 2. 



3 THE STAR FORMATION RATE 
EFFICIENCY IN HIGH-Z CLUMPS 

Equation ((9)1 shows that the most important parame- 
ter in determining the survival of giant clumps is how 
quickly, normalized to their free-fall times, they turn 
themselves into stars. Only if they do so with a very 
high rate efficiency, £9,-^2 ^ 10, do we expect signif- 
icant gas expulsion. In the local universe, one mea- 
sures eff by determining the star formation rate of an 
object or a population of objects, and comparing this 
to the objects' gas mass divided by the free-fall time 
computed for their density. Th i s pro cedure was first 
applied by IZuckerman fc Evans! (|l974l ) to the popula- 
tion of giant molecular clouds in the Milky Way, and 
has subsequently been exte nded to other obj ects by 
iKrumholz fc TanI (|2007l ) and lEvans et al.1 (|2009l ). These 
measurements give eft,_2 ~ 1 over a very wide range 
of star-forming environments, from small star clusters 
in the Milky Way to entire starburst galaxies. The re- 
sults are subject to considerable uncertainty, but values 
eff,-2 > 10 are strongly excluded by the data. However, 
we lack comparable data for star formation at high red- 
shifts. In this section, we therefore turn to the question 
of the likely value of eg in high-z giant clumps. 



3.1 Estimates of eg from Observations of 
Giant Clumps 

Unfortunately, we cannot easily apply the direct mea- 
surement procedure used to determine eg in the local 
universe to the high-z clumps, because, while we can 
evaluate star-formation rates using Ha luminosities, we 
are limited in our ability to measure the corresponding 

gas properties. 

For example, lElmegreen et al.l (|2009bl ) use stellar 
population synthesis to estimate masses and ages for 
the stellar populations seen in giant clumps in high-z 
galaxies, and they then compare the ages r to the clump 
dynamical times, defined as tdyn = {Gp)^^^^ ~ 0.5tff, 
where p is taken to be the stellar mass density. They 
find typical values r/tdyn ~ 1—10, corresponding to 
r/tff ~ 2 — 20, with a factor of ~ 10 scatter. It is 
tempting to identify r with t^cp and simply estimate 
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Eft ~ tfi/r ~ 0.05 — 0.5, but this is likely to be a sig- 
nifican t overestimat e . Bas ed on dynamical mass esti- 
mates, iGenzel et al.l (|2008l ) estimate that gas comprises 
10-30% of the total mass within the disc radius. Esti- 
mates based on direct CO measure ments indicate gas 
fractions of 45 - 6 0% at 2 ~ 2 jPaddi et all |2009| : 
iTacconi etall I2OI0I '). It is likely to be an even larger 
fraction of the mass within the dense, rapidly star- 
forming giant clumps, which are self-gravitating and 
lac k a dark matter component. If the stars measured 
bv lElmegreen et al.l comprise only a small fraction /, of 
the total mass in clump, then eg will be reduced by a 
factor of roughly /i'^ relative to the previous estimate - 
one power of /* to account for the gas that has not yet 
formed stars, and another factor of /"'^ because the free- 
fall time will be shorter than the value [Elmegreen et al.l 
estimate based on the stars alone. To give a sense of the 
ossib le magnitude of the error, note that iMurrav et al] 
20091 ) estimate /. = 0.2 for a giant clump in BX 482 
(and we argue in Section 14.21 that /, is probably even 
smaller) , and this value of /, would be sufficient to lower 
the estimate of eff by a factor of 1 0, to 5x 10~^— 5x 10~^ . 
A secondary worry is that, as lElmegreen et al] point 
out, since the clumps are selected using rest-frame blue 
light there is a strong bias against selecting older, red- 
der clumps, causing an underestimate of r. In general, r 
is expected to be an underestimate of tdcp because the 
former refers only to the stars that have already formed. 

Given this problem, many observers have at- 
tempted to esti mate gas masses via the "inverse" 
iKennicut^ (| 19981 ) law, namely by measuring the SFR 
surface density and then assuming that the gas surface 
density has the value required for the object to obey 
the Kennicutt law (e.g. lGenzel et al.ll2006l . l2008l 'l. Unfor- 
tunately this proce dure does n ot yield an independent 
estimate of eg . The iKenniciittl relation is 

S. = AHllt (10) 

with A ~ 1.4 Mq yr"^ kpc"^. Thus if the scale height 
of the galaxy is h, then the midplane gas density is p = 
S//i, the free-fall time is ts = ^37r/!,/(32GE), and we 
have 



es,^2 = 100 



-0.1 
-1 ) 



(11) 



where ho = ft/kpc and An = A/(1.4: Mq yr^^ 
kpc^^). The true value of eff is almost certainly a bit 
smaller than this, since star-forming clouds have densi- 
ties higher than the mean midplane density, and thus 
smaller free-fall times. Nonetheless, this calculation il- 
lustrates a crucial point: to the extent that galactic 
scale heights do not have a very large range of vari- 
ation (and the dependence is only to the 0.5 power), 
the statement that £ff.- 2 ~ 1 is roughly equivalen t to 
the Kennicutt law. (See Krumholz fc McKeell2005l and 



ISchave fc Dalla Vecchiall2008l for a more detailed discus- 
sion of the relationship between volumetric and areal 
star formation laws.) Thus, no measurement of the gas 
surface density that assumes the Kennicutt law a priori 



can produce a value of efr,_2 significantly different from 
this. Any measurement of eff,_2 that did yield a signif- 
icantly larger value would necessarily place the galaxy 
well off the Kennicutt relation. 



3.2 Estimates of eff from Observations of 
Clump Host Galaxies 

Given the difficulties of estimating eg in giant clumps 
directly, we instead turn to indirect inferences, based 
on the more robust measurements of the overall disc 
properties. We first note that the observed correlation 
between total gas mass and total star-fo r mation rate 
in high-2 galaxies (e. g. ICarilli et al. I l2005l :[ Greve et al.l 
l2005l : lGao et al.ll200iD is consistent with these galaxies 
havin g the same value of eg a s local star-forming sys- 
tems ( Krumhol z fc Thompson! l2007l : iNaravanan et all 
[2008a .b). Bothwell et all (|2009l ) have claimed to detect 
a deviation from a universal star formation law in spa- 
tially resolved_obseryations of three 2 ~ 2 systems. How- 
ever, |Bot^^irer^l]obt£^iii this result only because they 
choose a non-standard conversion factor between CO lu- 
minosity and mass, which leads them to conclude that 
the gas fraction in these systems is only ~ 5 — 10%, much 
lower than in typical star forming galaxies at 2; ~ 2, 
and on the low side e ven for local star-forming galaxies. 
iTacconi et all (|2010l '). using a standard conversion fac- 
tor, conclude instead that these galaxies have molecular 
gas fractions of 35 — 45%. This is consistent with the 
star formation law and the value of eft at 2; ~ 2 being 
the same as in the local universe. 

We can also approach the problem more theoreti- 
cally. The high-2 clumps discs with which we are con- 
cerned have SFR M* ~ lOOMgyr 
Md ~ 10" Mo, and disc radius _Rd 
express the relevant quantities in Equation ([9]) as a func- 
tion of the clump mass Af, surface density E, and SFR 
M,. The free-fall time (which is \/5/2 times the crossing 
time R/V), expressed in 10® yr, is 



^, baryonic mass 
- 10 kpc. We first 



iff, 6 



Then 



eff, -2 = 



3-K 

32Gp 



= 16.7A/g"^/^i?i''^ = 12.3Mg^'^T,Zl^'^ -{12) 



M, 
M/tti 



(13) 



where M^io is the SFR in a clump in 10 Mq yr"'^. 

Now let us express the clump quantities in terms of 
the disc mass, surface density, and total SFR. The disc 
surface density is 



Ed,-i = 0.663Md,iii?^- 



(14) 



where Rd,io is the disc radius measured in units of 10 
kpc. The virialized clumps can be assumed to have col- 
lapsed by at least a factor 2 in radius, so the clump 
surface density is s = 4s4 times larger than the disc's. 



E_i =2.65s4Md,iiJ?^jo. 



(15) 



© 2010 RAS, MNRAS OOO.nifTOl 



6 Mark R. Krumholz and Avishai Dekel 



with S4 ^ 1. 

Assume that a fraction a (~ 0.2) of the disc mass 
is in A*' (~ 10) identical clumps a nd a fraction /3 (~ 
0.5) of the SFR is in the clumps (|Dekel et al.l l2009bl : 
ICeverino et~al]|2009l '). Then 



M.,10 = 0.5/3o.5iVlo ^^.,100 : 

and 



M, 



9 



2Q0.2iVio'Afd,ll . 



(16) 



(17) 



Substituting the last three expressions in eq. p3p we 
obtain in terms of the disc quantities 



eff, 



-2 = 1.75s~'^*l3Q.5a-'^'''N-"''Ma*,iooMl-li''RX^^ .(18) 



-1/4, 



'3/2p3/2 



Note the very weak dependence on A''. We see that 
the most straightforward observational estimates for the 
massive clumpy di scs at z ~ 2 (e.g. the BzK galaxies, 
iGenzel et aLllboOST l yield efi,_2 1.75 and e ~ 0.1, i.e. 
no significant expulsion. 

This estimate is based on the observed properties of 
galaxies over a relatively narrow range of galaxy masses 
and redshifts, but using simple theoretical arguments 
we can deduce how the results are likely to scale to 
other galaxies for which we currently lack direct obser- 
vations. For a self-gravitating disc, the circular velocity 
is roughly ¥^200 - ^d,ii/-Rd,io, so eg -2 oc Md*V^^ , 
namely 



eff,-2 oc 



Afvir 



(19) 



A constant efi,_2 (~ 1) is consistent with the SFR be- 
ing a constant fraction of the baryon accretion rate, be- 
cause the la tter is roug hly proportional t o halo mass 
l|Neistein et a l. 2006, ; .Birnboim et all |2007| '1 : 



(20) 



The fact that the SFR follows the accretion rate is a 
natural result of the fact that the SFR is proportional 
to the mass of the available gas (|Bouche et al.l l2009l : 
iDutton et all |2009| ). and is consistent with the find- 
ing from simulation s when compared to observed SFR 
(|Dekel et aLlbOOgal ). We learn that e is only weakly de- 
pendent on M. 

The redshift dependence at a given halo mass, using 
cx A/v(l + zfl'^ and M oc (l-f zf , is 



eff,-2 oc s 



-3/4 



(1 + ^) 



1/2 



(21) 



The system can adjust the SFR to match the rate of 
gas supply by accretion with efi,_2 ~ 1 at all times by 
slight variations in the contraction factor s. At higher 
redshift, the clumps should contract a bit further and 
form stars at a somewhat higher surface density. 



3.3 Sub-Resolution Clumping 

Our discussion of clump survival in the preceding sec- 
tions is based on the assumption that giant clumps rep- 
resent single star-forming molecular clouds, although we 



of course expect them to possess significant substruc- 
ture, as do local molecular clouds. These substructures 
are unresolved by current observations and by simula- 
tions. Here we discuss how their presence affects our 
conclusions. 

First note that, for this purpose, we do not care 
whether any sub-clumps within the giant clumps them- 
selves survive star formation feedback and form bound 
stellar clusters. To see why, consider an extreme case in 
which all the sub-clumps within the giant clump expel 
most of their gas, and thus do not leave behind bound 
remnants. This is what we might expect to happen if 
all the sub-clumps had surface densities similar to that 
of their parent giant clump, but had masses well below 
the ~ 10^ — 10^ Mq minimum survival mass that we 
computed in § [51 In this case the sub-clumps would all 
form stars, expel their gas, and disperse, but both the 
stars and the expelled gas would still remain trapped 
within the much larger gravitational potential well of 
the giant clump. They could escape from this potential 
well only if the giant clump as a whole were disrupted 
by gas expulsion, which we have already shown in § [2] 
will happen only if eg is much larger than the expected 
value. Thus the end result of this scenario would be a 
bound giant star cluster without any bound sub-clusters 
inside it. 

At the opposite extreme, suppose that all the sub- 
clumps were to remain bound and undergo negligible 
gas expulsion. We might expect this scenario if the 
sub-clumps all had surface densities much higher than 
that of their parent giant clump. In this case the sub- 
clumps would convert most of their mass to stars, form- 
ing bound clusters. All the bound clusters would irra- 
diate the remaining mass in the giant clump, imparting 
momentum to it. If the stars imparted enough momen- 
tum, this gas would be expelled. Assuming most of the 
mass were in the inter-clump medium, as is the case for 
local molecular clouds, this expulsion would unbind the 
giant clump, producing many small individually bound 
clusters that are not bound to one another. Conversely, 
if the imparted momentum were not sufficient to un- 
bind the giant clump, as we expect, the result would be 
a giant star cluster consisting of many smaller bound 
clusters, all gravitationally bound to one another. 

In either extreme scenario, whether or not sub- 
clumps survive does not make any difference to whether 
a giant clump as a whole survives. This is dictated solely 
by the expulsion fraction from the giant clump. How- 
ever, sub-clumping still could make a difference for giant 
clump survival by raising the value of eff. In this case 
the sub-clumps would still have eff,_2 ~ 1, but the gi- 
ant clump would have ^^,-2 3> 1 because it would have 
the same star formation rate but a much lower mean 
density, and thus a longer free-fall time. 

To see whether this is likely to happen, we note 
that the turbulent motions within a giant clump are 
likely to break it up into smaller sub-clumps, much a lo- 
cal molecular clouds are broken up into clump, filamen- 
tary structures by turbulence. In such a configuration. 
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a majority of the mass is at a density higher than the 
volumetric mean density p that we have computed, and 
would therefore have a shorter free-fall time and a higher 
star formation rate. Quantitatively, we have computed 
the star formation rate as M* = egAf/fff (p), where M 
is the total mass of the giant clump, p is its volume- 
averaged density, and tg is the free-fall time computed 
at that density. However, if most of the mass is at a 
density p > p, the appropriate mass might be the mass 
M(> p) above that higher density, and the appropri- 
ate timescale might be ta (p) computed for that density. 
While one might worry that this cou ld be a significant 
effect, iKrumholz fc ThompsonI (|2007l ') point out that it 
is in reality quite small. Turbulent systems generally 
have lognormal density distributions. For such a dis- 
tribution, the fraction of the cloud mass with density 
greater than p is given by 



M{> p) _ 1 
M ^2 



1 + crf 



-21n3; -I- (7^ 

23/2(7^ 



(22) 



where x — p/p, p is the volumetric mean density, and 
(Tp is the dispersion of the density distribution. This 
is related to the Mach number M of the turbulence 
by CTp ~ l n(l -I- 'yM'^), where 7 is a constant of or- 
der u nity (|Padoan fc Nordlundl |2002| : iFederrath et all 
I2OO8I ). For (Tp = 2.5 — 3, the range of values expected for 
the Mach numbers found in giant clumps, the quantity 
M(> p)/ts[p) varies by only a factor of a few over a 
range of densities p/p « 10"^ - 10^. Thus even if most 
of the mass is at a density vastly larger than the mean 
density we have used, as long as the mass distribution 
follows the lognormal form expected for supersonic tur- 
bulence, the star formation rate will not be modified sig- 
nificantly from our estimate using the volume-averaged 
density. 



4 DISCUSSION 

4.1 Turbulence and Energy Balance in Giant 
Clumps 

Our finding that the fraction of gas ejected from giant 
clumps depends critically on their dimensionless star- 
formation rate efficiency ea naturally leads to the ques- 
tion of how this quantity is set, and whether the phys- 
ical processes responsible for setting efi,_2 ~ 1 in the 
local universe m ight determine a diffe r ent va lue in high- 
redshift clumps. IKrumholz fc McKed (120051 ) show that, 
as long as the gas in a molecular cloud is supersonically 
turbulent with a velocity dispersion comparable to the 
cloud's virial velocity, as is observed to be the case in 
all molecular clouds in the local universe, efF,-2 ~ 1 is 
the inevitable consequence. In contrast, in the absence 
of supersonic turbulence, simulations find that clouds 
undergo a rapid global collapse in which they convert 
all their mass into s tars in roughly a dynamical time, 
i.e. eff,_2 ~ 100 (e.g. iNakamura fc Li|[2007l : IWang et al.l 



|2009| )FI Thus, a value of eff,-2 ~ 1 may be expected in 
high-2 giant clumps only if they maintain the level of 
turbulence required to avoid rapid, global collapse. 

Whether the turbulence can actually be maintained 
is somewhat less clear. Simulations show that supersonic 
turb ulence decays in r oughly one cloud-cro s sing time 
(e.g. IStone et al.ill998l : iMac Low et al.lll998l : iMac Lowl 
I1999I ). so global collapse can be avoided only if this 
energy is replaced on a comparable timescale. In lo- 
cal, lo w-mass star-formin g clouds (< 10 "* Mp)), observa- 



tions JOuillen et al.ll2005l l. simulations (INakamura fc L: 
2007 : 1 Wang et all 120091 ). and analytic theorv (|Matzner 



20071 ) all suggest that protostellar outflows can supply 



the necessary energy. In local giant molecular clouds 
with masses of ~ 10^ — 10^ Mq, H 11 regions driven 
by the pressure of photoionized gas a re likely to be 
able to supply the necessary energy (|Matzneij |2002| : 
IKrumholz et al.ll2006l ). However, neither of these mecha- 
nisms are effective for clumps with Afg ~ 1 and E_i ~ 1, 
because they do not provide enough momentum input 
and because the y are overwhelm ed by radiation pressure 

(see Figure 2 of lFall et al1l2010l). 

Supernova feedback (|Dekel fc Si ii3 ll986l ) does not 
appear to be a likely candidate to drive the turbulence 
either. Supernovae do not provide e nough power to driv e 
the observed level of turbule nce (iDekel et al.ll2009al). 
and analytic calculations ( Harper-Clark fc Murravl 
I2OO9I : IKrumholz fc Matznerl |2009|). numerical simula- 
tions of isolated disk galaxies ( Tasker fc BrvanI l2008l : 



Ijoung et aLll2009l ). and numerical simulations of galax- 
ies in cosmological context (|Ceverino fc Klvpinl l2009) 
all indicate that supernova-heated gas is likely to escape 
through low-density holes in the molecular gas without 

driving much turbulence. 

Contrary to this conclusion, iLehnert et al.l (|2009l ) 
use the observed correlation between Ha surface bright- 
ness and linewidth in z ~ 2 galaxies to argue that super- 
nova feedback is responsible for driving the turbulence, 
based in part on simulations bv lDib et al.l (|2006l ). who 
obtain a scaling relation between velocity dispersion and 
supernova rate in numerical simulations. However, the 
efficiency with which supernova energy is couple d to the 
ISM is a free param eter in both lLehnert et aLr s analy- 
sis and in lDib et al.l 's simulations, and their results are 
consistent with the data only if it is ~ 25%, whereas in 
the dense environments found in high redshift galaxie s 
it is expected to be far lower (|Thompson et aLll2005l ). 
This conclusion is confirmed by the more recent simula- 
tions, which do not need to assume an efficiency because 
they have sufficient resolution to resolve the multiphase 
structure of the ISM. Finally, we note that the corre- 
lation be tween Hq surfa ce brightness and linewidth ob- 
served bv lLehnert et al.) has a more prosaic explanation: 



* T ills can be avoided if c louds are magnetically subcritical 
(e.g. INakamura &: Li|[20o3 ). but magnetic fields in the early 
universe are likely to be weaker than those in the local uni- 
verse, and even in the local univ erse clouds do not ap pear to 
be subcritical in typical clouds (ICrutcher et al.ll2009l) . 
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in a marginally stable galactic disk of constant circular 
velocity, the velocity dispersion is proportional to the 
gas surface density, since Q = na / {nGYl) = 1. Thus 
higher velocity dispersions correspond to higher surface 
densities, which in turn produce hig her star formation 
rates in accordance with the standard iKennicuttI (| 19981 ) 
relation. This naturally explains the observed correla- 
tion. 

Radiation pressure is another mechanism to con- 
sider. If radiation pressure is not able to drive mass out 
of the clump, as found above for es ~ 0.01, this sug- 
gests that it might not be able to drive turbulence to 
the required virial level either, since the virial and es- 
cape velocities only differ by a factor of ^/2. However, 
it is unclear whether this conclusion is warranted. Ra- 
diation pressure cannot drive material out of a clump 
not because stars do not accelerate material enough, but 
because they evolve off the main sequence before they 
are actually able to eject matter. As a result, they pro- 
duce expanding shells whose velocities greatly exceed 
the escape velocity. They simply fail to drive mass out 
because the clump because the driving sources turn off 
before the shells actually escape from the cluster. It is 
unclear if the expanding shells might provide enough 
energy to maintain the turbulence; this problem will re- 
quire further modeling. 

We are left with the possibility that the turbu- 
lence is driven by gravity. The driving source cannot 
be the collapse of the clump itself; although such a col- 
lapse does produce turbulence, it does so at the price 
of reducing the crossing time, raising the rate of en- 
ergy loss. Consequently, the collapse becomes a run- 
away process, and all t he gas quick l y conve rts to stars. 
However, as shown by iDekel et al.l (|2009al ). the gravi- 
tational migration of the clumps through the galactic 
disc does provide enough power to maintain the turbu- 
lence within them. The main uncertainty in this model 
is how much of that power will go into driving internal 
motions within the clump, rather than motions in the 
external galactic disc. This depends on how the clumps 
are torqued by one another and by the disc. However, 
there is suggestive evidence from simulations of giant 
clumps that this mechanism might be viable. The sim- 
ulations of clumpy galaxies that have been done to 
date (e.g. Bournaud et all 20071: Elmegreen et alllioOSl : 



lAgertz et al.ll2009l : ICeverino et all 120091 ') either include 
no feedback or only supernova feedback (which is inef- 
fective). If the giant clumps formed in these simulations 
did lose their turbulence and undergo global collapse, 
then, depending on the details of the simulation method, 
they would either convert all of their mass into stars, or 
all of the mass within them would collapse to the maxi- 
mum density allowed by the imposed numerical pressure 
floor. This collapse would happen on the crossing time 
scale of a clump, which is much less than the time re- 
quired for the clumps to migrate to the galactic centre. 
However, such collapses are not observed in the simula- 
tions. This strongly suggests that, even in the absence 



of feedback, gravitational power is sufficient to maintain 
the turbulence. 

We conclude that the generation of turbulence in 
the giant clumps is an important open issue, to be 
addressed by further studies including simulations of 
higher resolution. 



4.2 Comparison to Previous Work 

Our conclusion that giant clumps are not likely to be 
disrupted by feedbac k is in contrast with the findings of 



iMurrav et al.l (|2009|) for the giant clump in the galaxy 
Q2346-BX 482 (|Genzel et al.ll2008l '). and this difference 



merits discussion. Ba sed on the Ha luminosity of the 
clump, iMurrav et al.l estimate a total bolometric lumi- 
nosity of L = 4 X 10^^ Lq, which corresponds to a star- 
formation rate of 34 Mq yr~^ in the old stars limit, 
or a stellar mass of 2.6 x 10* Mq in the young stars 
limit, which iMurrav et al.l assume. If these stars are in- 
deed young, then the star-formation rate could be higher 
than 34 Mq yr~^, but it could not be any lower. 

While the estimate of the star-formation rate is rel- 
atively straightforward^^ the gas mass is much 
less so. iMurrav et al] (|2009l ') take it to be 10^ M© based 
on an order-of-magnitude estimate for the Toomre mass 
in the galaxy. This choice is crucial to their result. The 
clump radius is r = 925 pc, so if we adopt this radius, 
the mean density-free fall time is 15 Myr, so es,-2 — 50. 
Using this value in equation ([HI, together with the cor- 
responding surface density and mass E_i = 0.8 and 
Mg = 1, tells us that we are in the y oung stars limi t 
and that £ = 0.06 - fully consistent with IMurrav et all 's 
conclusion that only a relatively small fraction of the 
gas mass turns into star s, and that th is is sufficient to 
expel the remaining gas. [Murray et al.l derive a slightly 
higher value f ~ 1/3 because their criterion for ejection 
amounts to adopting Qcrit ~ 10. 

However, this conclusion depends crucially on hay- 
ing a low estimate of the gas mass, and a correspond- 
ingly high estimate for eff,_2. Indeed, if the value of 
es,-2 — 50 were accurate, this clump would have the 
highest star formation rate efficiency of any known sys- 
tem. It is therefore useful to consider alternative meth- 
ods for estimating the mass. If we were to adopt the 
inverse-Kennicutt method, the observed star-formation 
rate per unit area E* ~ 13 Mq yr~^ kpc~^, together 
with equation pO|l . gives a gas surface density of E_i ~ 

5 (or 2300 Mq pc~^). The corresponding gas mass is 

6 X 10® Mq, and recomputing eg for this mass gives 
eff__2 — 3, consistent with the point we made earlier 
that the Kennicutt Law is in practice equivalent to hay- 
ing eff,-2 ~ 1. With this mass we would predict £ ~ 0.9, 
i.e. essentially n o gas expulsion . One would expect sim- 
ilar results from [Murray et all 's models, because the ef- 
fect of this mass increase would be to increase the grav- 
itational force by a factor of 40 while leaving t he radia- 
tive force unchanged. Thus while lMurrav et al.l 's models 
suggest that the clump in BX 482 has stopped forming 
stars and all the remaining gas has just been expelled. 
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according to our estimate this clump can be only part 
of the way through its life and may continue to form 

starS; 

iMurrav et ai] (|2009l ') also suggest another method 
of estimating the gas mass. Based on the Ha luminosity, 
if one assumes that the clump is filled with uniform- 
density gas, then ionization balance requires that this 
gas have a density of hydrogen nuclei 

"^ = ^4^3^(5) (5^)' (23) 

where X is the hydrogen mass fraction and this ex- 
pressio n assumes tha t He is singly-ionized. If, fol- 
lowing I Murray et al.l (|2009l ). we adopt the young 
stars limit, then the ionizing luminosity is S = 

6.3 X W^^jMJMp^ ) = 1.6 X lO'^^ photons s"^ 

(|Murrav fc RahmanI I2OO9I ') Fl Combining this with the 
case B recombination coefficient q'^' = 3.46 x 10"^"^ 
cm^ s~^ and the Solar hydrogen mass fraction X = 0.71 
(mean mass per H nucleus of lAmu) gives nn ~ 21 
cm~ , corresponding to a mass of 2.4 X 10® Mq and 
eff,-2 = 14. Plugging this mass, surface density, and 
value of eff,_-2 into equation ((S)) gives £ ~ 0.2, i.e. the 
star fraction is more than t hree times what we would 
obtain using [Murray et al.l 's mass of M = 1 x lO'' 
Mq. Adopting Ocrit = 10 to s hift our fiduc ial pa- 
rameters closer to those used in [Murray et al.l would 
give £ ~ 0.8, no significant gas expulsion. We empha- 
size that these calculations are lower limits on the gas 
mass and upper limits on the fraction of mass ejected, 
because this mass estimate includes only ionized gas. 
However, models of both classical gas pressure-driven 
H 11 regions and ones driven by radiation pr essure 
(|Krumholz fc MatzneJ l2009l : [Murray et all \200^ sug- 
gest that the ionized gas mass in the H 11 region in- 
terior is significantly smaller than the mass of neutral 
gas swept up in the shell around it. Including this mass 
would lower efi^_2 and increase £ even further. 

In summ a ry, ou r conclusions differ from those of 
[Murray et all (|2009l ) not because of any difference in 
the physics of radiation feedback, but because they have 
used an estimated gas mass that produces an extraor- 
dinarily high value of eff. 



5 CONCLUSIONS 

Our main result in this paper, summarized in Equation 
((5]) and Figure [T] is that the survival or disruption of 
giant star-forming clumps in high-2 galaxies depends 
critically on the rate at which they turn into stars. We 
find that as long as the high-redshift clumps convert 
their gas mass into stars at a rate of one to a few per- 
cent of the mass per free-fall time, ~ 0.01, as is ob- 
served in low-redshift star-forming systems from small 

^ The ionizing luminosity is somewhat lower in the old stars 
limit: 5 = 2.5 x 10^^{M,/Mq yr^^) = 8.5 X 10^* photons 
s-i. 



galactic clusters to ultraluminous infrared galaxies, the 
clumps retain most of their gas and turn it into stars. 
As a result, they remain bound as they migrate into the 
galactic centre on timescales of ~ 2 disc orbital times. 
A significant fraction of the clump gas could be ejected 
on a free-fall timescale before turning into stars only if 
clumps can convert > 10% of their gas mass into stars 
in a free-fall time, forming stars much faster than any 
other star-forming system known. We argue the current 
high-2: data is consistent with the standard SFR rate 
efficiencies at the level of one to a few percent, with 
no significant evidence for a change in the star forma- 
tion process in high-z star-forming galaxies. Neverthe- 
less, this is clearly an interesting issue to explore with 
more direct observational estimates of the gas mass in 
these galaxies. 

It is possible to check our theoretical arguments for 
clump survival with a number of possible observations. 
First, clumps can be disrupted by radiation pressure 
only if a majority of their gas is expelled, and the result- 
ing massive radiation-driven outflows from clumps may 
be observable as systematic blueshifts at the clump loca- 
tions. A preliminary search for such a phenomenon have 
yielded a null result (K. Shapiro, private communica- 
tion, 2009), but further investigations of the kinematics 
in and around the clumps are worthwhile. It is possible 
that sufficient extinction could hide the blueshifted sig- 
nature, but we note t hat, if the gas s urface densities are 
relatively low as, e.g.. [Murray et all (|2009l ) propose, the 
extinction is relatively mild. 

Second, stellar populations in high-redshift clumps 
can provide independent observational tests that could 
help distinguish between the two scenarios of clump 
survival or disruption. One such test involves the age 
spread of stars in actively star-forming clumps. If ra- 
diative feedback disrupts the clumps after one or a few 
free-fall times, the spread of stellar ages in each clump 
should not exceed ~ 50 Myr. If clumps survive, on the 
other hand, then the age spreads may reach > 100 Myr, 
with a high and roughly constant SFR during the life- 
time of the clump. There is p reliminary observational 
evidence in favor of the latte r ([Elmegreen et al1[2009bl : 
[Forster Schreiber et al.l [20091 ) , but this should be ex- 
plored further. In this case, migration toward the bulge 
will produce an age gradient, so that clumps closer to 
the bulge have systematically larger age spreads than 
those further out in the disc. 

A third test is associated with the properties of 
massive star clusters in the disc that are not actively 
forming stars. Clusters with masses > lO'^ Mq can only 
be made in giant clumps, not in the rest of the disc, 
where molecular clouds are smaller. If clumps undergo 
rapid disruption by feedback, most of the clusters are 
likel y to dissolve as is the case for local star formation 
(e.g. [Fall et"ai][2009l '). but the few clusters that may sur- 
vive will remain in the disc near their formation radii. 
None will migrate into the bulge, sin ce the migration 
time varies as jPekel et al ] [2009b[ ). and the masses 
of the clusters are much smaller than those of their par- 



© 2010 RAS, MNRAS OOO.ITlfTOl 



10 Mark R. Krumholz and Avishai Dekel 



ent gas clumps. In contrast, if gas clumps survive feed- 
back and convert most of their mass to stars, the re- 
sulting massive objects will rapidly migrate toward the 
galactic center. Th ey will lose ~ 50% of t heir mass due 
to tidal stripping (|Bournaud et al.ll2007l ). possibly in- 
cluding some massive sub-clusters, but they will deliver 
the rest of their mass to the bulge. Thus, if giant clumps 
do not survive, we expect to see a few massive clusters 
in the disk and none in the bulge. If, on the other hand, 
our model is correct, then there should be comparable 
masses of disk and bulge clusters. Those bulge clusters 
that survive today may correspond to the metal-rich 
globular cl uster population seen in present-day galac- 
tic bulges (|Brodie fc StradeJbOOd IShapiro et allboiOl : 
Romanowsky et al., in preparation). 

Even though we have concluded that the high-z 
giant clumps are expected to survive radiative stellar 
feedback, we do emphasize the important role that this 
process is likely to play in galaxy formation. For exam- 
ple, it can provide some of the pressure support needed 
in these giant clumps, it can drive non-negligible winds 
out of them, and it is likely to disrupt the less mas- 
sive clumps where stars form at later redshifts. In many 
circumstances, this mode of feedback is expected to be 
more important than supernova feedback, as it pushes 
away the cold dense gas while the latter mostly affects 
the hot dilute gas. We thus highlight the need for in- 
coroprating radiative stellar feedback in hydrodynam- 
ical simulations as well as in semi-analytic models of 
galaxy evolution. 
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